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Abstract 

We present a systematic cosmological study of a universe in which the visible sector is 
coupled, albeit very weakly, to a hidden sector comprised of its own set of particles and 
interactions. Assuming that dark matter (DM) resides in the hidden sector and is charged 
under a stabilizing symmetry shared by both sectors, we determine all possible origins of 
weak-scale DM allowed within this broad framework. We show that DM can arise only 
through a handful of mechanisms, lending particular focus to Freeze-Out and Decay and 
Freeze-In, as well as their variations involving late time re-annihilations of DM and DM 
particle anti-particle asymmetries. Much like standard Freeze-Out, where the abundance 
of DM depends only on the annihilation cross-section of the DM particle, these mechanisms 
depend only on a very small subset of physical parameters, many of which may be mea- 
sured directly at the LHC. In particular, we show that each DM production mechanism 
is associated with a distinctive window in lifetimes and cross-sections for particles which 
may be produced in the near future. We evaluate prospects for employing the LHC to 
definitively reconstruct the origin of DM in a companion paper. 



1 Introduction 



As our understanding of theoretical high-energy physics has evolved, top-down considerations 
have motivated the exploration of "parallel sectors" comprised of their own particles and inter- 
actions but "hidden" from us (the visible sector) due to the weakness of the couplings connecting 
hidden and visible sector particles. In addition to providing new avenues for model-building, 
this broad framework opens up a range of exciting possibilities for the origin of Dark Matter 
(DM), which is the focus of this work. Understanding the origin of DM and its interactions 
within this framework is very important because experimental observations have only measured 
the gravitational effects of DM, leaving a large number of logical possibilities. 

Assuming that the visible sector and other possible hidden sectors are initially in a state 
of thermal equilibrium, what are the possible production mechanisms for DM? If DM shares 
sizeable interactions with visible sector particles, then thermal equilibrium will be efficiently 
maintained until Freeze-Out (FO) renders a thermal relic abundance of DM via the standard 
WIMP paradigm [I]. Alternatively, it may be that DM couples extremely weakly to the visible 
sector and to itself, as is the case for so-called superWIMPs [2j [3] and FIMPs [U |5] . A third 
and final possibility is that DM is very weakly coupled to the visible sector, but has substantial 
couplings to a hidden sector to which it is thermally equilibrated. In general, this hidden sector 
will contain its own set of particles and interactions and will have a temperature different from 
that of the visible sectoiQ. The purpose of the present work is to systematically identify and 
characterize all possible origins of DM which might arise in this enormous class of theories. 

We will assume throughout that DM is stable due to a symmetry shared by the visible and 
hidden sectors. Moreover, let us denote the lightest visible and hidden sector particles charged 
under this symmetry by X and X', which have masses m and m' taken to be broadly of order 
the weak scale such that m > m! . By definition, X' is the DM particle. We also assume the 
existence of a weak coupling which bridges the visible and hidden sector and mediates the decay 

X^X' + ..., (1) 

where the ellipses denote what are typically visible decay products. 

Remarkably, the cosmological evolution of this setup is entirely fixed by only a handful of 
parameters. This is analogous to standard single sector FO, where the DM abundance is solely 
determined by the DM annihilation cross-section. Here we find that DM relic abundance is fixed 
by following set of parameters in general: 

{m,m',{ov),(ovy,€,T,e}, (2) 

If the temperature of the two sectors are the same, the sectors have equilibrated implying that there is only 
one sector. 
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where (av) and (crv)' are the thermally averaged annihilation cross-sections for X and X', 
respectively, £ is the ratio of the visible and hidden sector temperatures, r is the lifetime of 
X, and e is a measure of the CP-phase in X decays. In particular cases, the relic abundance 
depends on only a subset of the above parameters, as will be shown below. 

We have evolved the cosmological history of the visible and hidden sectors over the parameter 
space defined in Eq. (J2J) in order to systematically identify all possible origins of hidden sector 
DM. Of course, the simplest possibility is that DM undergoes hidden sector Freeze-Out (FO'), 
yielding a thermal relic abundance. This has been considered in many hidden sector models, 
and was studied systematically in [6J. On the other hand, the remaining possibilities for the 
origin of DM fall into two very broad categories: 

• Freeze-Out and Decay (FO&D). X undergoes FO and then decays out of equilibrium, 
yielding an abundance of X'. As we will see later, the final abundance of X' goes as 

~ ml 

n oc -j—. (3) 

m{av) 

• Freeze-In (FI). X decays while still in thermal equilibrium with the visible sector, yielding 
an abundance of X'. As we will see later, the final abundance of X' goes as 

n oc ^L. (4) 

Within the categories of FO&D and FI exist a number of distinct variations. For example, if 
FO&D or FI happen to produce an abundance of X' particles exceeding a particular critical 
value, then the X' particles will promptly undergo an era of "re-annihilation." During this 
time the X' particles will efficiently annihilate within a Hubble time despite the fact that X' 
is no longer thermally equilibrated with the hidden sector. Because the final DM abundance 
changes accordingly, we refer to this mechanism of DM production as FO&D r and FI r . Another 
variation arises if X decays are CP-violating, in which case FO&D and FI may produce an 
abundance of DM endowed with a particle anti-particle asymmetry. Such an effect is possible 
because although the visible and hidden sectors are separately in thermal equilibrium, they are 
not in equilibrium with each other. We denote these asymmetric modes of DM production by 
Asymmetric Freeze-Out and Decay (FO&D a ) and Asymmetric Freeze-In (FI a ). Note that these 
mechanisms are entirely distinct from the framework of Asymmetric DM [TJ, in which the DM 
particle anti-particle asymmetry is inherited from an already existent baryon asymmetry. 

Crucially, as seen in Eqs. |3] and 0J each of these DM production mechanisms maps to a 
rather distinctive window in the parameter space spanned by r and (erf) — and where all other 
parameters, m, m' , (crv)', £, and e, are scanned over an inclusive range of values. This is 
remarkable because r and (av) can, in principle, be measured at the LHC — after all, they are 
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Figure 1: Hidden sector DM can originate via a handful of production mechanisms, each corresponding to a 
distinctive window in the r — (av) plane. Aside from hidden sector FO', these mechanisms are {FO&D, FO&D r , 
FO&D a , FI, FI r , FI a }, denoted by {blue, green, purple, red, orange, yellow}. Each point corresponds to ilh 2 = 
0.11, where we have scanned over a very inclusive parameter space defined by 10~ 5 < (av) / (av)o, (av)' / (av)o < 
10 5 , 10 -3 < £ < 10 _1 , 10~ 8 < e < 10~ 3 , where (av) = 3 x 10~ 26 cm 3 /s. In the left panel, the masses have been 
scanned over a broad region 10 GeV < m < 1 TeV and 1/20 < m'/m < 1/2, while in the right panel, the masses 
have been fixed to a narrow region m = 100 GeV and 1/4 < m'/m < 1/3. 

attributes of X, which is a visible sector field! For example, see the left panel of Figured! where 
each point corresponds to Qh 2 = 0.11, and each color denotes the dominant mechanism of DM 
production at that particular point in parameter space. Even though all parameters but r and 
(av) have been scanned over a generous range, one sees that FO&D corresponds to a narrow 
band in (av ) while FI corresponds to a narrow band in r. 

On the other hand, it is also very likely that m, and perhaps even m', might be measured at 
colliders, for instance if the visible decay products of X — > X' + . . . can be used to kinematically 
reconstruct the event. In the event that m and vn! are indeed measured, the boundaries between 
DM production mechanisms in the r — (av ) plane become even more distinct, as shown in the 
right panel of Figure [TJ 

Because each production mechanism lies in a distinctive region in the r — (erf) plane, we are 
left with the tantalizing possiblity that the origin of DM might be successfully reconstructed at 
the LHC even in this much broader framework compared to that of standard single sector FO. 
The purpose of the present work, however, is to establish a comprehensive understanding of the 
structures depicted in Figure [fl leaving a more detailed collider study to a companion paper [S] . 

The outline of this paper is as follows. In Section |5] we present a broad overview of two- 
sector cosmology. We begin with an analysis of our setup in a decoupled limit in which the 
visible and hidden sectors couple only through gravitational interactions. We then introduce 
portal interactions, and present a detailed discussion of the FO&D and FI mechanisms of DM 
production. Afterwards, we go on to discuss the effects of "re-annihilation", followed by an 
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analysis of the thermal properties of the coupled two-sector system. In Section [3j we present a 
series of cosmological phase diagrams depicting the dominant production mechanisms for DM 
as a function of the parameter space. We go on to discuss how the boundaries in these phase 
diagrams change with various parameters. In Section H] we present a discussion of DM production 
from particle anti-particle asymmetries, and we conclude in Section [5j 

2 Overview of Two-Sector Cosmology 

Our setup is comprised of a visible and hidden sector, each with sizeable self-interactions which 
serve to maintain thermal equilibrium in each sector at temperatures T and T", respectively. 
We assume that these sectors couple to one other only through portal interactions which are 
extremely feeble, so these temperatures are not equal, i.e. T ^ T' . To begin, we limit the 
present discussion as well as that of Sections 12.11 and 12.21 to the case in which the visible and 
hidden sectors are entirely decoupled but for gravitational effects. In Section I2T51 and onwards we 
introduce portal interactions connecting the visible and hidden sectors and study the significant 
impact of these couplings on the cosmological history. 

Throughout, we assume that the visible and hidden sectors enjoy a symmetry, discrete or 
continuous, that keeps DM cosmologically stable. The lightest visible sector particle charged 
under this stabilizing symmetry is denoted by X, and likewise in the hidden sector, X', which 
we take to be lighter than X. In the limit in which the visible and hidden sectors are decoupled, 
X and X' are, of course, simultaneously stable. However, as portal interactions are switched on, 
X becomes unstable and decays with a width T into particles which ultimately yield an X' in 
the final state. We will study the cosmological evolution of the number densities n and n! of X 
and X', which obey the coupled Boltzmann equations 



where n eq and n' eq are the thermal equilibrium abundances and it is understood that here T is 
thermally averaged and we work in a regime in which the effects of the corresponding inverse 
decays are negligible. Here we take the thermally averaged annihilation cross-sections, (av ) and 
(crv)', to be independent of temperature. A primary aim of this paper is to study the most general 
cosmological evolution which follows from these equations, subject only to the requirement that 
(av) and (av )' are large enough that both X and X' undergo freeze-out. Note that we take the 
masses of X and X', m and m', both to be broadly of order the weak scale. 

The relative size of T and T" can have a drastic impact on the cosmological history. To see 



— n + 3Hn 
at 




(5) 




(6) 
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why this is so, let us define the ratio of temperatures to be 

V 

i = y- (7) 

If we assume the standard picture of slow-roll inflation, then the inflaton can, in principle, couple 
with different strengths to the visible and hidden sectors. As a consequence, the decay of the 
inflaton reheats each sector to a different temperature, corresponding to an initial condition for £ 
given by £i n f = T( nf /Ti n f, the ratio of temperatures in each sector immediately after the decay of 
the inflaton, which we take to be less than 1. Naively £ = £i n f for all time. However, interactions 
between the sectors can change £ from £ in f. For example, scatterings between the sectors which 
are generically dominated in the UV, can increase the high temperature value of £ to £uv> which 
is taken to be a free parameter. In addition, there can be IR contributions to £ as well. These 
contributions are discussed in more detail in Section 12.71 

Even ignoring interactions between the two sectors, £ actually varies as a function of temper- 
ature due to the separate conservation of the co- moving entropies, S = g*sT 3 and S' = g'^ s T' 3 , in 
each sector. Specifically, this implies that £ varies as a function of temperature to the extent to 
which the numbers of degrees of freedom in the visible and hidden sectors vary with temperature: 



where g*s{T) and g^ s (T) are the number of relativistic degrees of freedom in the visible and 
hidden sectors, respectively, when the visible sector is at a temperature T. A change in g*s(T) 
or g'^ s {T) by an order of magnitude only affects £ at the level of a factor of two; hence, when 
comparing DM production from processes at different temperatures, this effect may be justifiably 
ignored. 

In general, the energy density in the hidden sector affects the expansion rate of the universe 
during BBN, which places an important, albeit weak constraint on £. In particular, any hidden 
sector particles which are relativistic at BBN contribute an effective number of extra neutrino 
species 



The present bound from experiment is AN U < 1.4 [9], which is surprisingly mild: that is, for 
S*C?bbn) = 100, tlris is satisfied by taking £(T B bn) = 1/3. Furthermore, according to Eq. ([8]), at 
higher temperatures £ can be close to unity even if g[ > 100 [6]. 

2.1 Visible Sector Freeze-Out (FO) 

In the early universe, visible sector particles reside in a thermal bath at temperature T with 
abundances fixed accordingly by equilibrium thermodynamics. As T drops below the mass of X, 
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AN U = % ^(Tbbn) £(Tbbn) 4 . 



(9) 
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m, the number density of X particles, n, remains in thermal equilibrium and undergoes the usual 
Boltzmann suppression. The X particles undergo FO as the rate of annihilations, n(o~v), drops 
below the expansion rate H, which occurs at a temperature T F q — m/x-po- The parameter x-po 
depends only logarithmically on (crv), and for roughly weak scale cross-sections x F o ~ 20 — 25. 
Defining the yield, Y = n/s where s is the entropy of the visible sector, we obtain the familiar 
expression for the yield at FO, 

Y ~ A fl ^E i i (1Q) 

FO 27rV2 g* s M Fl {av)T FO ' 1 ) 

In the decoupled limit, X is stable and will account for the totality of DM in the universe if 
mYpo ~ 4 x 10~ 10 GeV. This corresponds to a critical cross-section of (crv)o - 3x 10~ 26 cm 3 /s. 

Because we have thus far assumed that the visible and hidden sectors only interact gravita- 
tionally, the only effect of the hidden sector on FO in the visible sector is through its contribution 
to the energy density of the universe. However, this effect is tiny and can be accounted for in 
Eq. (flDl by increasing g* by a factor 



i ,i 

3 



43 " \g*{T BW ) \9l(T F o) 



where we have ignored the difference between g* and g*s- For T FO < M w this increases Ypo by 
at most 17% for AiVj, = 1. Thus the BBN constraint implies that the standard relation between 
the DM abundance and the DM annihilation cross-section is preserved to a good accuracy. 

2.2 Hidden Sector Freeze-Out (FO') 

As the temperature T' of the hidden sector falls below m', the number density of X' particles, 
n', tracks the equilibrium distribution and becomes exponentially suppressed. Ultimately, the 
X' particles undergo FO' once the hidden sector drops to a temperature Tp , = m'/x' FO ,, just 
as X undergoes FO at Tpo = m/x F o- Here the prime on T' indicates that the temperature is 
that measured in the hidden sector, while the prime on FO' indicates that this temperature is 
being evaluated at the time of FO', not FO. Thus, the ratio of visible sector temperatures at 
FO' compared to FO is 

^ = ~ (12) 

where we took x' FO , ~ xfo- Consider the parametric scaling of the above expression. Reducing 
m! jm tends to shift FO' to a later time than FO, while reducing £ FO ' tends to do precisely the 
opposite. While in principle either ordering is possible, we focus on a scenario in which FO' 
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Figure 2: X and X' yields as a function of m/T for a particular choice of parameters. The visible and hidden 
sectors separately undergo FO and FO', respectively. At high temperatures Y'/Y ~ £y V = 10~ 6 . Since ml is not 
much less than m, a cool hidden sector, £ <C 1, implies that FO' occurs before FO. A cool hidden sector tends to 
make Y^q, < Ypo, but this is more than compensated by having (crv)' <C (crv). Using the quoted parameters in 
Eq. CEO) gives Y^ Q , ~ 200 Y fo . 

occurs before FO, as this leads to a richer set of cosmological histories once portal interactions 
between the sectors are included. 

First, consider the case {av)' = {crv). Because n' at FO' is fixed by H at T F( y and T FO ' > Tpo, 
there is naively a greater number of X' particles yielded by FO' than X particles yielded by FO. 
However, because FO' occurs earlier, there is also a commensurately greater amount of entropy 
dilution, so the total X' yield actually turns out to be less than the X yield. To see this, let us 
define Y' = n' / s to be the X' yield normalized to the visible sector entropy, which will be useful 
for comparing with the X yield. We find that for arbitrary {av)' / {av) 



Hence, the yield of X' from FO' is sub dominant to the yield of X from FO as long as the hidden 
sector is sufficiently cool or if its annihilations are sufficiently strong. See Figure [2] for a plot 
of the evolution of X and X' abundances as a function of x — m/T, for a choice of parameters 
making the contribution from FO' greater than that from FO. 




(13) 
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2.3 The Portal 



Until now, we have not considered the effect of direct, albeit tiny couplings which might directly 
connect the visible and hidden sectors. Consider a portal operator O which connects X and X', 
thereby mediating the decay 

X ->• X' + ..., (14) 

where the ellipses denote what is typically visible SM particles. For the moment, let us ignore 
the particulars of O and attempt to characterize the gross features of the cosmological history 
as a function of the X lifetime, r = 1/T. As the lifetime is taken from cosmological scales to 
microscopic scales, the cosmology typically transitions through four broadly defined scenarios^: 

• Multi-Component Dark Matter. X is so long lived that it is stable over cosmological 
time scales. Thus X and X' comprise the DM of the universe. 

• Freeze-Out and Decay. X decays late, after leaving thermal equilibrium, yielding a 
contribution to the X' abundance. 

• Freeze-In. X decays fast enough that it produces a substantial X' abundance from decays 
occurring while X is still in thermal equilibrium. 

• Thermalized at Weak Scale. X decays so quickly that the visible and hidden sectors 
are actually in thermal equilibrium at the weak scale. From the point of view of cosmology, 
the visible and hidden sectors are a single sector. 

While the first category is certainly a logical possibility, it has been well explored in the literature 
and is hard to test experimentally since the DM abundance depends on (cry )', so we will ignore 
it. Moreover, we will not consider the last category because we are specifically interested in 
cosmological scenarios in which the visible and hidden sectors are not thermally equilibriated at 
the weak scale. Thus, our discussion will center on the FO&D and FI phases of the two-sector 
cosmology. 

2.4 Freeze-Out and Decay (FO&D) 

In the presence of the portal operator, O, X is no longer stable. Thus, after X undergoes FO, 
it eventually decays into X' particles; we call this DM production mechanism "Freeze-Out and 
Decay" (FO&D). The resulting X' may form the dominant contribution to the final yield of X', 
as illustrated in Figured Assuming the X decay process, X — > X' + . . ., produces exactly one 
X' for each X, we find 

^FO&D - *FO- (15) 
2 This is only a rough sketch; a more precise understanding of the various possibilities is given in Section [H 
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Figure 3: X and X' yields as a function of m/T for DM production dominated by FO&D. The X particles 
undergo FO and then later decay, yielding an X' abundance that is larger than that which arises from FO'. The 
parameters £, (av),m and m! are the same as in Figure [2 but (ov)' is increased giving Yp , > Ypo- For r = 1 
second, X decays are occurring at the MeV era. 

Consequently, the energy density produced by FO&D is suppressed relative to that of conven- 
tional FO by a factor of m'/m. If FO&D accounts for the total DM abundance in the universe, 
then this implies that (av) = (m! /m)(av)o, where recall that (crv)o is the annihilation cross- 
section needed to account for the measured DM abundance in standard single sector FO. This 
dilution factor is useful in theories in which FO normally produces an overabundance of DM, 
for instance as occurs in supersymmetric theories if the LSP is a bino. 

Mechanisms similar to FO&D have been discussed extensively in the literature for a small 
subset of candidates for X and X' and operators O. In particular, there is a large body of 
work [21 [3] concerning the so-called superWIMP scenario in which X is effectively a bino or 
right-handed slepton and X' is the gravitino. Axinos [10] and goldstini [Tlj [12] have also been 
studied as alternative choices for X' . 



As the X decay rate is increased, at a certain point a new DM production mechanism, "Freeze- 
In" (FI), begins to dominate. Here X' particles arise from decays of X particles which are still 
in thermal equilibrium. As long as the X' have already undergone FO', the X' produced by FI 
can comprise the dominant source of DM, as shown in Figure HJ 

At any temperature T > m the production of X' by FI generates a yield which goes schemat- 



2.5 Freeze-In (FI) 
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Figure 4: X and X' yields as a function of m/T for DM production dominated by FI. While X is relativistic and 
in thermal equilibrium, the small fraction of X that decay yield an important contribution to the X' abundance. 
As soon as X' undergoes FO', the FI mechanism begins to effectively populate an X' abundance that grows 
until T drops below m, when the X abundance becomes exponentially suppressed. Once the age of the universe 
reaches r, the relic X particles from FO all decay but, for the parameter choice for this figure, the increase in Y' 
from this FO&D process is sub-dominant to the FI contribution. 



ically as 



Y Fl {T) oc Tt oc 



pi 



J^2 



(16) 



where t is the total time that X is relativistic. A key aspect of FI by decays is that it is 
IR dominated by low temperature dynamics; this is true independent of the dimensionality of 
the connector operator which mediates the decay. FI can also occur by two-to-two scattering 
via a marginal coupling (this is also IR dominated). However it turns out to be numerically 
subdominant compared to that from decays and inverse decays [3], so this will not be discussed 
from now on for simplicity. As T becomes larger, FI plays an important role in increasing £(T) 
as the temperature drops, as we will discuss in Section 12.71 Here we focus on the X' produced 
after FO'. The FI yield from X decay is dominated by contributions from T ~ m and is the 
same as computed in [I] for FI from inverse decays. The precise formula for the FI yield is 

TMpi 



V 



Cfi(^fo' 



Cfi(^fo' 



135 15 g x 



(17) 



Ki(x)x 3 dx 



FO ^ ' 



Fl 



1.64- 



9x 



2vr 5 V ^g*s^/9* Jx FO , 9*sV9* 

where gx is the number of degrees of freedom of X and x FO ' = m/T FO i. As shown in Eq. [12], in 
the limit in which the hidden sector is much cooler than the visible sector, X' freezes out very 
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Figure 5: X and X' yields as a function oim/T for DM production dominated by FI followed by re-annihilation. 
Re-annihilation occurs because the FI yield exceeds Y^ lit at some temperature where 1q SE > Y^ lit . Once this 
happens, Y' tracks Vqse until it dips below Y^ lit , which occurs as Y drops rapidly at FO. The value of Y' at this 
crossing point is the final yield of X' particles. 



early so Xpo' — ► 0. Finally, note that if FI accounts for the total DM abundance today, then for 
weak scale masses this implies a range of lifetimes given by r ~ 10 -4 s — 10 _1 s. If the decay 
of X is mediated by a marginal operator with the dimensionless coefficient A, then this range 
of lifetimes corresponds to A = 10~ 12 — 10~ n . For decays mediated by a higher dimensional 
portal interaction, this range applies to A = (m/M„) d ~ 4 where d and M* are the dimension and 
scale of the higher dimension operator. Interestingly, for d = 5 operators, this corresponds to 
M* ~ 10 13 — 10 15 GeV, which is roughly of order the GUT scale. 

2.6 Re- Annihilation 

For simplicity, we have ignored the effects of X' annihilation on FO&D and FI. Naively, this is a 
justifiable omission, since both FO&D and FI occur only after the hidden sector has undergone 
FO'. However, even after FO', the X' abundance arising from non-equilibrium production may 
be so large that the X' annihilation rate grows to exceeds the expansion rate, initiating a new era 
of X' annihilation that we dub re-annihilation. Re-annihilation and the resulting X' abundance 
can be understood through a study of the Boltzmann equation, written in terms of yield variables 
and x = m/T, 

x d Y' - Y ' 2 + TY (19) 

dx ^crit H 
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where 

C - ^ (20) 

The first term on the right-hand side of Eq. [19] corresponds to X' annihilation; since we are 
interested in times well after FO', Y' 3> Y£ and we ignore inverse annihilations. The second 
term is effectively a source term for X' production, corresponding to the decays of X to X' . At 
T ~ m this is the source term which drives FI, while for T w y/TMpi this is the source term 
which drives FO&D. However, the analyses of FI and FO&D in the previous sections ignored 
the annihilation term. 

The destruction and production of X' occur faster than the Hubble rate if the first and 
second terms on the right-hand side of Eq. f fl9|) are larger than Y', respectively, that is if 



Y'>Yr t & ir >Y '- (21) 



In this case the Y' abundance rapidly evolves to a Quasi-Static Equilibrium (QSE) in which the 
production of X' particles is counter-balanced against depletion from the annihilation process. 
This causes the two terms on the right-hand side of Eq. ffl9|) to cancel, so that Y' becomes equal 

to *qse> where 

/2 TY , TY 

YqSE = -JfYcrit = J^y~ s - (22) 

Setting Y' = Yq SE in the first equation of ( 12T1) . one discovers that QSE is possible only during 
eras having 

^QSE > Krif (23) 

As X' undergo FO', the depletion of X' will stop once Y' drops to Yq SE , provided Eq. fl23|) 
is satisfied, as shown in Figures |5] and |6] which were produced by numerically solving the exact 
Boltzmann equations. Subsequently Y' tracks Yq SE until Eq. ( 12 3 p is violated. This always 
eventually happens because Y^ lit grows linearly with m/T and 1q SE drops as Y is reduced by FO 
or X decay. When QSE ends, i.e. when Yq SE = Y" c ' rit , the re-annihilation rate drops below the 
expansion rate so that QSE is lost and the X' yield becomes constant. The final X' abundance 
is fixed by the value of Y' at this point. We will denote FO&D and FI which are subsequently 
followed by re- annihilation by FO&D r and FI r . Note that re- annihilation did not occur in the 
examples shown in Figures |2]|4] because Eq. (|23|) was violated at FO' and all subsequent eras. 

Figures and show yield plots for cases where DM is dominated by FI r and FO&D r , 
respectively. The dashed lines indicate Yq SE and Y^ rit . In both plots, one sees that, once 
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Figure 6: X and X' yields as a function of m/T for DM production dominated by FO&D followed by re- 
annihilation. Y' tracks Yq SE until it drops below Y^ lit , which occurs when Y drops sharply as X decay. 



inverse annihilations of X' can be neglected, QSE occurs during eras with Yq SE > Y^ it , with 
Y' accurately tracking Yq SE . Eventually, YA SE falls below Y^ Tit and QSE ends. The final DM 
abundance reflects the value of Yq SE at the end of the QSE era, which we now study. 

To analytically compute ^po&D r an d Y Flr , we must compute the temperature at the end of 
QSE, T r , which is by definition the solution to the equation 



^qse(^t) 



^crit(-^r)- 



(24) 



According to Eq. (|22j) . Yq SE oc y/Y, so the solution to this equation depends on the expression 
for Y during the era under consideration. Going from T to x variables, then for FO&D r and FI r 
we must solve the transcendental equations 

V4 



-T 1 / 2 p-^FO&Dr — I 

x FO&D r e ~ 2 

5/ 2 -x FI _ 



187T 2 g* 

V27T 7 / 2 g* 



g M^r(av)' ' 



FO&D r 
FI r 



(25) 



Hence, the final yields for FO&D and FI followed by re-annihilation are given by 



V 

FO&D r 

Y' 



^crit(^FI r ), 



(26) 
(27) 



This result is very similar in structure to the yield obtained from standard FO', namely Y FO , = 
y^. it (T FO ')- For s wave annihilation Y" c ' rit oc 1/T, so the re-annihilation yields obey the simple 
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relation 



Tfo&d.^fo&d, - ^Fi r ^pi r - T FO/ Y FO , {21 
To be concrete, this implies that the DM yield for FO&D r and FI r are given by the formulas 



^FO&D r 2?r y 2 g ^ ^ ( ay y Tpo&Dr 

FI " 2n\ 2 g* s M Pl (av)>T Flr ' 1 ; 

where T F i r = m/x F \ r , and T F o&D r = ^Decay/i/^FO&D,. so that the exponential in Eq. (125]) goes 
simply as e~ x . Here T Decay is the temperature at which X decays and xpi r and x F o&d t are given 
by the solutions of Eq. (125]) . 

As we have seen, the DM yield from FO&D and FI can differ substantially from FO&D r 
and FI r . The condition for avoiding re-annihilation effects is Ypj < Y" c ' rit (T = m) for FI and 
Y FO < F c ; it (T = v/TMpT) for FO&D. 



2.7 Sector Equilibration 

Thus far we have ignored the effects of the connector operator O on the thermal properties of 
the visible and hidden sectors. Specifically, there is the danger that O couples the sectors so 
strongly that they actually come into thermal equilibrium. This scenario corresponds to the case 
where £ = T'/T ~ 1 near the weak scale. In general, £ is temperature dependent and receives 
contributions from UV and IR-sensitive physics, 

m = £uv + £i 4 R (T). (30) 

The UV contribution to the hidden sector temperature, £uv ; arises from two sources, so £{j V = 
£inf + Cr- If f ne inflaton couples directly to the hidden sector, then £ in f is generated by an initial 
heating of the hidden sector from inflaton decays. This contribution is independent of the portal 
interactions, and was discussed in earlier sections. On the other hand, £r results from scattering 
processes mediated by O occurring at the reheat temperature, T R . These processes are active if 
O is a higher-dimension operator. In this case, O contributes a 2-to-2 scattering cross-section, 
(crv) R , which produces an X' yield of 

M Pl T R (av) R . (31) 

The T R dependence of (crv) R depends on the dimensionality of O. If X' is inert, Y R can easily 
overdose the universe unless T R is sufficiently small. For instance, in the case of gravitino LSP 
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this is the origin of the well-known bounds on T R from overclosure [13]. On the other hand, 
our assumption is that X' possesses self-interactions, so X' particles produced by scattering at 
reheating will be efficiently thermalized via the X' annihilation until the onset of FO'. Only 
after FO' can an abundance X' particles be produced via FI. Consequently, in the presence of 
X' annihilations, the FI abundance from the higher dimension operator O is given by Eq. (3TJ 
only with Tr replaced by Tpo'- Because T FO ' is not exceedingly far from the weak scale, this 
UV dominated FI contribution from 2-to-2 scattering will in general be subdominant to the IR 
dominated FI contribution from decays discussed in Section 12.51 

Since the X particles are produced with energy T R , the hidden sector is reheated by the 
visible sector to an energy density given by T' R 4 ~ Y R T R . Thus, the ratio of visible and hidden 
sector temperatures is given by 

^ R ~ (M Pl T R (av) R y/\ (32) 

in the case where there is 2-to-2 scattering processed mediated by higher dimension operators. 

Finally, let us consider the IR contribution to £, which essentially arises from FI. For decays 
of X at temperature T, FI produces a yield Yp T given by f|T6|) . The produced X' particles have 
an energy distribution characteristic of temperature T. We assume that the interactions of the 
hidden sector are sufficient to rapidly thermalize the energy of these X' into distributions of all 
the hidden sector particles at temperature T". As long as the hidden sector remains sufficiently 
cool, then FI will be mediated by decays and inverse decays can be ignored. From this point of 
view, FI leaks energy and entropy out of the visible sector into the hidden sector. This effect is 
especially important when the lifetime of X is short, since the strength of this energy leakage is 
proportional to T. The contribution to £ir from FI is 

Zm(T) = A^f- (T>m), (33) 

with an analytic estimate giving A m 135v5 gx/(ir 5 y/thg'*)- As T drops below m the FI process 
gets exponentially switched off, so 

6R(T)~£ m (m) (T<m). (34) 

As T increases we reach a critical point where £ir(T ~ m) = 1; the FI process is now so strong 
that the two sectors are equilibrated when T ~ m. Since m > m' the X particles no longer 
undergo FO, and instead track their equilibrium abundance by decaying. Consequently, there 
is single thermalized sector and DM results from the single process of FO'. For £uv ^ 1> the 
critical lifetime that leads to this equilibration scales as m~ 2 

100GeV\V 100 



r min ~ 1(T 13 s ( ) I — ) , (35) 

m J \9*{T ~m) g x ' 
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where the numerical coefficient has been extracted from our numerical results. Since the contri- 
bution to the heating of the hidden sector from FI is model-independent, we include its effects 
throughout our analysis. On the other hand, the UV-sensitive contributions are very model- 
dependent, so we take £uv to be a free parameter which is small. 

2.8 Summary of Results 

Let us now summarize the results of the entire section. We have shown that in our setup 
the present day abundance of DM particles can originate only from a handful of cosmological 
production mechanisms. In the simplest case, the hidden sector undergoes FO', yielding a 
thermal relic abundance of X' particles. Alternatively, via the FO&D mechanism, an abundance 
of X particles can FO in the visible sector, and then decay very late into X' particles. The FI 
mechanism functions so that X particles, while still in thermal equilibrium, provide an abundance 
of X' particles through decays. Lastly, if the X' yields from FO&D and FI exceed a certain 
critical yield, Y^ rit , then the hidden sector enters an era of re-annihilation. The final abundances 
for FO&D r and FI r are controlled by the temperature at the end of this re-annihilation era. The 
analytic expressions for the DM yield are 

^fo&d = Cfo -rj—, — — oc — - — - (36) 
Mpi{av) T FO m{av) 

Mpi(av)' T FO ' m'(avy 

^Fo&D r — Cfo irj~/ — \7t^ ^ ( 3 ^) 

Yk = C FO M ) ~ oc — ^— (39) 
Mpi{avyTpi t m{avy 

= C Fl ^ oc — 2 (40) 

in the cases where the various mechanisms dominate the contribution to the total yield. Here 
we have defined the dimensionless constants Cpn = ~n~ \ / n and Cft = 1.64 — a -=, and the 
various temperatures are defined in Table [U The dimensionless values u x" defined in the second 
column of this table are determined in each case by solving a transcendental equation of the 
general form: 

x n e~ x = f{m,m',(av),(av)',Z,T) (41) 

for some rational number n and where / is a some function of the arguments. Here we have 
taken an approximate solution in which the effect of x n is neglected. 
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Table 1: Expressions for the various temperatures relevant for each DM production mechanism. The "x" 
quantities employed in the first column are given approximate expressions in the second column. The quantity 
£fo' is computed in (|42l) below. 



Note that only FO' depends on £; in particular it depends on the value of £ at T FC y which 
is denoted in the Tabled] as £ F o'- From the analysis in Section [277] one notes that the quantity 
£ F( y = £(T FO ') has different forms depending on whether T F( y is greater or smaller than m. £ F( y 
can be computed as: 



;fo' 



_|_ ATMpi^l/4 



ATM Pl x 



2 m'- 



1/2 



1 + 



4$ v m ' 4 

A 2 r 2 m 2 



PI ^FO 



1/2' 



1/2 



T FO ' < m 
T FO ' > m 



(42) 



where A is as defined after Eq. ( J33l . 



3 Cosmological Phase Diagrams 

A primary aim of this paper is to identify and characterize all possible mechanisms of DM 
production which can arise within our general two-sector framework. To this end, we have 
simulated the cosmological history of this system over a broad range of values for the relevant 
parameters: 

{m, m', (cry), (av)', £uv, r}, (43) 

where £tjv is the UV initial condition for £ which receives contributions from the decay of the 
inflaton as well as scattering processes from higher- dimensional operators described in Eq. (IBTJIl R 
As noted earlier, it is quite remarkable that the cosmology is determined solely by just a handful 
of quantities. 

3 Here we also take £uv to include effects from additional sources of entropy dumping into either sector before 
the weak era, so that £uv is effectively the weak scale value of £, modulo the contribution from X decays in the 
IR. 



17 



In this section we present a series of "cosmological phase diagrams" depicting the regions in 
parameter space in which each mode of DM production, i.e. FO&D, FI, etc., accounts for the 
dominant contribution to the present day DM abundance. For example, Figure [7] is a cosmolog- 
ical phase diagram in the r — ^4- plane for particular values of {m, m', (crv), £uv}, as explained 
in the caption. Each colored shaded region corresponds to a particular DM production mode 
which is dominant in that region. The boundaries of each phase have been computed analyti- 
cally. The solid black contours correspond to the present day DM relic abundance, calculated 
numerically using the coupled Boltzmann equations in Eqs. (J5J). Let us consider some of the 
features of Figure [7] in detail. 

First, note that the lifetime r is constrained both from above and below as shown by the 
red vertical lines in Figure [7J The upper limit is easy to understand as it originates from the 
requirement the decay products of X do not ruin the successful predictions of BBN [13] . While 
the precise constraint depends on the nature of X and its decay products, here we use 

W ^ 100 s, (44) 

as depicted by the red line at large r in Figure [7J The lower limit on r comes from demanding 
that the two sectors are not in thermal equilibrium with each other at the weak scale. In Section 
12.71 we found that for £y V <C 1, the two sectors becomes thermally equilibrated if the lifetime is 
shorter than 

r ~ 10-" s ( 100GeV ) 2 ( ^ 1 (45) 

For r < 10~ 13 s, the two sectors are, from a cosmological perspective, a single sector. As 
a consequence, there is no distinction between FO and FO' — the visible and hidden sectors 
together maintain thermal equilibrium until X' undergoes single sector FO'. This is depicted 
by the brown region in the left in Figure [7J The yield Yp Q , in this region is independent of r 
because £po' = 1 in the expression for Y^ Q , in Table [TJ Note that the brown region extends to r 
larger than r min , as will be discussed soon. 

Next, let us discuss the salient features of Figure [7J as we move from the largest to the 
smallest allowed values for r. Furthermore, for a given value of r, the dominant DM production 
mechanism changes as a function of 7^4-- For example, for the choice of parameters in Figure 
[7J one finds that FO&D and FO' dominate over FI for 1CT 1 s < r < 100 s. This is because 
the lifetime is too long (the coupling is too weak) for the FI mechanism to set in, as shown in 
Section 12.51 For values of r in the above range with ^4- > 10 _1 , corresponding to the blue 
shaded region, FO&D is the dominant mechanism. Indeed, from the shape of the contours one 
can see that in this region the relic abundance is independent of ^4- and r as dictated by the 
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Figure 7: Cosmological phase diagram showing regions in the versus t plane where different mechanisms 

contributing to the relic abundance dominate. The values of the other relevant parameters are chosen as: 
£ uv = 0.01, m = 100 GeV, m! = 50GeV, (av) = (av) = 3x 10" 26 cm 3 /s. Contours of Vth 2 , computed 
from a full numerical analysis for Y' (numerical solution of Eq. ([5]) and Eq. ©), are shown. Regions in which 
the various mechanisms dominate are shown in different colors. These regions are computed analytically and 
overlaid on the numerical plot. The agreement is quite good. 

expression for Ipo&D i n Tabled] For ^4- < 10~ 2 , FO' starts to dominate over FO&D since the 
annihilation cross-section of X' becomes sufficiently small. This is shown by the brown region 
in the right in the figure. For such long lifetimes, one finds that £ FO ' ~ £uv from (!30l) . implying 
that Yp Q , is essentially independent of r (see the expression for Yi Q , in Table [p. Thus, the 
contour lines in this region are roughly horizontal. The boundary between the FO&D and FO' 
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regions for this range of r is given by : 



fo&d fo' =► - &C^7 



oc constant (46) 



(aw) 

which in the ^4 — r plane corresponds to ^4- ~ constant, since £ FO ' ~ £uv- 

For 10 -8 s < r < 10~ 4 — 10 _1 s, FI and FI r begin to dominate over FO&D and FO' because 
the portal interactions between sectors is growing stronger. In particular, FI is dominant in the 
narrow band shown by the red shaded region in Figure [71 For FI to be the dominant mechanism, 
the yield from FI must be larger than that from FO' but smaller than that from FI r ; hence the 
narrow band in which FI is the dominant mechanism. The boundary curves between the FI 
region and the FO' & FL. regions are given by: 

Y > - Y> =± EMM C ™ T 



(<TV) 

r2 



oc t (47) 



TM^(av)'m' _ C FO 



(avY 

(<TV) 



OC T (45 



Both curves imply a linear relation between ^4- and r, only with different coefficient^ 4 ). Also, 
note from Table [T] that Fpj does not depend on ^4-, hence the contour lines for the relic 
abundance are almost vertical. They are not completely vertical since it turns out that FI never 
fully dominates the contribution to the relic abundance. 

For 10~ 8 s < r < 1CT 4 — 10 _1 s with t^4- < 1, the FI yield becomes so large that it exceeds 
the critical yield at Tpi r , implying that X' starts to re-annihilate and the yield is given by Yp X . 
This is shown by the orange region in the figure. From Table (TJ one sees that contours for Y Flr are 
essentially horizontal since the dependence on r only arises from the logarithm. As one increases 



above roughly unity, FO&D starts to dominate over FL, as shown by the blue region, since 



(avy 

(<TV) 

the annihilation cross-section of X' becomes large enough that the yield Y FIr becomes smaller 
than Y FO& , D . The boundary curve between these two regions is defined by: 

(av)' x Fh 



V — V 

-'FO&D _ -'FI,, 



av) x FO 

<TV)'1 



(av)' 

oc Xpj oc In 



(av) 



(av) r 



(49) 



4r 



The boundary shown in Eq. (|48|) depends on £po' as well, but in this region £po' ~ £uv is a constant. 
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showing that is essentially constant up to a logarithmic dependence. 



If 



Ml 

(crv) 



(av) 

is further increased, FO&D is eventually superseded by FO&D r as shown by the green 



region. In this case the X' yield from the freezeout and decay of X is so large that it becomes 
larger than the critical yield; so X' starts re-annihilating with a yield given by i^o&D ■ From 
Tabled], the contour plots for Yp Q ^ Dr follow a simple power law oc y/r). The boundary 

curve between the FO&D and FO&D r regions is given by: 



Y' 



V' 

FO&D 



(avY 
(av) 
(avY 
(av) 



1 



V^FO&Dr 



m 



^FO&D r ^fo a/TMpi 



'OC 



(50) 



showing that the boundary also follows the same power law as the DM yield contours. 

For few x 10~ 13 s < r < 1CT 8 s, FL is again the dominant mechanism as shown by the orange 
region. It dominates over FO&D r in particular, since the yield is inversely proportional to the 
relevant temperatures for FO&D r and FL, which are T F o&D r and T FIr respectively, and T FOr is 
larger for small lifetimes (large T) compared to T F i r from Table [TJ The boundary curve between 
these two regions is given by: 
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I FI r 
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(51) 



(av) 



and r 



which is roughly r oc constant, although more precisely there is a dependence on 
from the logarithms. 

Finally, for r min < r < few x r min , FO' again dominates for the entire range of as is 
shown by the brown region in Figure [71 To understand this, it has to be compared to the next 
relevant mechanism, namely FI r . Since the lifetime r is short in this region, £fo'> which is given 
by the expression in the top line in Eq. f|42|) in this case, becomes large and close to unity. Hence, 



in this region, £ F( y > tllEEI^^ or equivalently T F( y < T F i r < m from Table CD Since (av) has been 



fixed in Figure [7J from Table [TJ and 

(av)' 



the contours for the DM yield in this region have a 



dependence ^v- oc r x / 4 . The boundary curve between the two regions is given by: 



*FO' 
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(52) 
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Figure 8: Sequence of cosmological phase diagrams with (crv) varying. Other parameters are as in Figure [7] 
(£ uv = 0.01, m — 100 GeV, ml = 50 GeV). The black contour corresponds to flh 2 = 0.11, the observed DM relic 
abundance. 



which is again approximately r oc constant, but more pecisely with an additional dependence 
on and r in the logarithm. 

((TV) ° 

3.1 Behavior of Phase Space Diagram 

In the cosmological phase space diagram depicted in Figure 7, the parameters {m, m', (crv), £tjv} 
were fixed. Let us now examine how each phase region changes as we vary these four parameters. 
Figures [HI El and [TU] are cosmological phase diagrams showing the variation of one of these 
parameters while keeping the remaining three fixed. The effects on the regions can be understood 
by referring to the formulae in Section [3] which define the boundaries between regions. In 
addition, we would like to examine how the mechanism of DM production is affected. To this end 
we include in these figures a single "critical DM abundance" contour (black line) corresponding 
to the Qh 2 = 0.11. 

Figure [H] illustrates the effect of varying (crv). The most dramatic change is in the contours of 
total yield. As can be seen from the analytic yield formulas in Section [31 when (crv) is decreased, 
must increase in order to maintain a constant yield, whether dominated by FO', FO&D, 
FO&D r , or FI r . In addition the boundaries of the FI region shift to smaller lifetimes. Thus the 
yield contours "rise" in the r — ^4- plane. For larger values of (erf), as in the left panel of 
Figure El the black line corresponding to Qh 2 = 0.11 can access the FI region but not the FO&D 
region. Once FI dominates it will give the correct DM abundance when: 

r ~ (3xlO- 15 s)xM P1 -A (53) 

m gj 

m'\ /100GeV\ /228.5 X 3/2 



(4xir2 (54) 
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Figure 9: Sequence of cosmological phase diagrams with £uv varying. Other parameters are as in Figure [7] 
(to = 100 GeV, to' = 50GeV, (av) = 3 x 10~ 26 cm 3 /s). The black contour corresponds to flh 2 = 0.11, the 
observed DM relic abundance. 



this corresponds to a lifetime of ~ 1CT 2 s in the left panel of Figure [H As (av) is decreased 
the critical DM contour rises and can start to access the FO&D region. For FO&D to give the 
correct DM abundance the following relation must hold: 

(av)m 4 x 10 10 2 x 1(T 25 cm 3 sec _1 . . 

■ ' - (55) 



m 1 M Ply /Sk jgl 

This is the case for the center panel of Figure [8] with g± = 228.5. Decreasing (av) even further 
leads to FO&D r domination, as shown in the right panel of Figure El 

Figure M illustrates the effect of varying £tjv- As is expected from Eq.(46) decreasing £uv 
results in a smaller large-r FO' region. This makes sense since decreasing £uv corresponds to an 
earlier FO' (for large r, £fo' ~ £uv) an d a smaller Y^ Ql which will thus dominate in a smaller 
region of parameter space. It is clear that the FI region expands as £uv is decreased. 

Figure [10] illustrates the effect of varying ml . From (|54|) it is clear that the lifetime that 
yields the observed DM abundance from FI is proportional to m! . In the left panel of Figure [T0| 
w! = 50 GeV and the FI lifetime is r ~ 10~ 2 s; decreasing w! it is possible to lower this lifetime. 
This is shown in the center and right panels of Figure [10J m! = 10 GeV requires r ~ 10 -3 s 
for FI to dominate. Also note that decreasing m! results in a larger FO' region at small r, as 
expected since this boundary is approximately the line £ = — . 

Finally we comment on the variation of m. The only major effect is the change in the lifetime 
at the thermalization bound (red line) which corresponds to £ (r, m) = 1. Thus as m is increased 
this lifetime is decreased. 
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Figure 10: Sequence of cosmological phase diagrams with ml varying. Other parameters are as in Figure [7] 
(£ uv = 0.01, to = 100 GeV, (av) = 3 x 10~ 26 cm 3 /s). The black contour corresponds to Slh 2 = 0.11, the observed 
DM relic abundance. 

4 Dark Matter Asymmetry 

Until now our discussion has been limited to cases in which the final abundance of the DM, X', is 
symmetric under the exchange of particles and anti-particles. As such we have implicitly assumed 
that either X' is its own anti-particle or that the mechanism of DM production dynamically 
produces equal numbers of X' and X'. In this section we determine the conditions necessary for 
DM to arise from a particle anti-particle asymmetry 

. n' — n' . . 

rf = 56 

s 

rather than from a symmetric yield, Y' = (n' + n')/s. An asymmetric mechanism of DM 
production requires that the hidden sector possess a global U(l) symmetry, Q'. A crucial 
question is then whether the portal interactions connecting each sector either completely break 
Q' , or whether they preserve some combination of Q', baryon number B, and lepton number 
L. We study two possible scenarios in which there exists a global U(l) symmetry 5* having the 
properties 

• Si = Q' rotates only hidden sector fields and is broken by the connector interactions. 

• S 2 = aB+f3L+ r yQ' , with real parameters a, (3 and 7, is exact, except possibly for anomalies. 

4.1 Asymmetric FI and Asymmetric FO&D 

At temperatures well above the weak scale we assume that there does not exist any asymmetry 
between the DM particles and anti-particles, so 77' = 0. Is it possible for a non-zero asymmetry 
to be generated by the thermal production mechanisms studied in the last section? At first 
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sight the answer is no, since each sector is separately in thermal equilibrium. While thermal 
equilibrium is lost right at the end of FO and FO', these processes involve total annihilation cross- 
sections (cry), (cry)' which are the same for anti-particles as for particles, so an asymmetry cannot 
be generated. However, the sectors are at different temperatures, so processes mediated by the 
connector interactions are not in thermal equilibrium. This offers new possibilities for generating 
a DM asymmetry, and indeed for baryogenesis. For the symmetry Si, the connector interactions 
may generate rf, but it is unrelated to the baryon and lepton asymmetries, T)l,b- On the other 
hand, symmetry S 2 is exact so any Q' charge that is generated must be compensated by baryon 
and lepton charges, so that the universe remains neutral under the charge of S% ■ Asymmetric FI 
that provides such a correlation between the baryon and DM densities is discussed in [14J, and 
explicit supersymmetric models are constructed. 

There are now four Boltzmann equations to be studied for the evolution of n, n, n' and n' . 
The equations for n and n are identical and take the form of Eq. The decay rate T is the 
total decay rate of X, and by CPT the total decay rate of X is also T. However, the equations 
for n' and n' become more complicated once we consider multiple decay modes for X, which is 
one of the necessary conditions to generate an asymmetry. We index each decay mode by i, with 
a partial width Tj corresponding to 

X ->■ piX' + piX' + (particles with Q' = 0) (57) 

where Pi and pi are integers denoting the number of X' and X' particles produced via the ith 
decay mode. The X' and X' may result from a cascade of decays in the hidden sector. In general, 
there can be modes containing particles other than X' and X 1 (say Y') with Q' 7^ 0; however we 
assume that the hidden sector interactions are large enough that Y' rapidly decays into X' and 
X', so that we only need to consider modes of the form in fl57j) above. If X is not equal to its own 
anti-particle then X has a set of decay modes to the corresponding anti-particles with partial 
widths Fj, and by CPT the total decay rates of X and X are equal, i.e. defining T = Tj and 
F = ^\ Fj, then r = f . On the other hand, if X is real it can decay both via Eq. f|5T|) and by 
the decay to the corresponding anti-particles, so that its total lifetime is T = Tj + £\ Fj. 
The Boltzmann equations for n' and n' are: 

—n' + 3Hn' = -{n'n' - n'^n'^) (cry)' + n^pjTj + n ^pjti (58) 

i i 
i i 

As with the analysis of the symmetric abundance Y' in the previous section, we omit scattering 
process contributions to FI. They are UV dominated both for Y' and n', and it is unclear how 
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they can be distinguished experimentally from a high scale initial condition on rj' '. Taking sums 
and differences of these two equations, the source terms for n! ± n' are Y2i(Pi ^Pi)( n Ti ± nTi)- 
The goal of this analysis is to understand how rf is generated from these equations starting from 
symmetric boundary conditions, i.e. when n = n. As such, this mechanism is markedly different 
from using these equations to transfer a pre-existing asymmetry in X to X', as has been studied 
in the context of ADM [7J. With n = n the source terms become Fn^^pi ±pj)(rj ± fj), where 
we introduced the branching ratios = I\/r and fj = fj/f . This last form is also the correct 
source term for the case that X is the same as its anti-particle. 

Using these source terms, the symmetric and asymmetric contributions to the DM yields are 
given for FO&D by 



Mpilpoiav) — ' MpiiFoiav) z — ' 

i i 

and for FI by 

*fi = Cfi—^- y2(Pi+Pi)( r i + f i) Vfi = Cfi^^- Vfe - Pi)(r< - f<). (61) 



The symmetric yields Y FO& . D and Y Fl were studied in great detail in Sections |2] and [3j Note 
that these yields do not take re- annihilation into account; this will be discussed soon. On the 
other hand, the asymmetric yields rj' FO& . D an d Vfii are new — we denote these DM production 
mechanisms by Asymmetric Freeze-Out and Decay (FO&D a ) and Asymmetric Freeze-In (FI a ). 
As seen from (IBDl) and ([ST]) , generating an asymmetry requires decay modes with pi ^ p, so that 
the final state has Q' ^ 0, as well as fj ^ 

According to Eqs. (1601) and flBTT) . the symmetric and asymmetric mechanisms obey the general 
relation, 

rf = eY\ (62) 

where 

e = >p (gj -pi)(rj -fj) 
~ (Pi+P^i^ + ^i)' 

Here e is a general measure of CP violation occurring in decays of X. If, for example, there are 
two relevant decay modes such that p\^ and p\^ are not large, then e is roughly given by r\ — f\. 

Next, let us determine the typical size of e for a simple case with two decay modes. We take 
the two decay modes of interest to be X — > X' + f\ and X — > Y' + f 2 , where fi are particles 
in the final state that have Q' — 0, and Y' and X' have different Q' charge. We introduce 
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dimensionless amplitudes A4 to describe these decays, defined by Tj = |Aj| 2 m/87r. Furthermore, 
as mentioned earlier we assume the existence of additional (hidden) interactions which allow for 
the rapid decay of Y' via Y' — > p 2 X' + p 2 X', with p 2 — p 2 7^ 1- In addition, a vertex allowing for 
the rescattering process Y' + f 2 — > X' + fi (with dimensionless amplitude A12) is also required. 
This allows the final state of process 2 above to rescatter into the final state of process 1 at one 
loop (and vice versa), which is necessary for successful asymmetry production. 

If the aforementioned amplitudes are too large, they will cause the two sectors to equilibrate 
at the era T m and destroy the viability of the FI and FO&D mechanisms. Applying the 
equilibration condition Eq. fl35l) demands that 



\A I < in" 6 / m &{T~m)/g x 

VToOGeVV 100 • (64) 

A non-zero value for e results from interference between tree and one loop contributions to the 

decays 

_ 1 lm(A 1 A* 2 A 12 ) (g5) 



16vr l^l 2 + \A 2 \ 2 ' 

In general the rescattering involves both visible and hidden sector particles so that, to avoid 
equilibration of the sectors at T ps m, A\ 2 must satisfy the same bound, ( |64|) . as Aj, giving 



e < 10- 8 sin <P J ^J 9 ^"^ 9 * (66) 
* V 100 GeV V 100 V ; 

where = arg^i^A^) . 

Is it possible to evade this bound? For theories with a global symmetry of type S2 such that 
a combination of B, L and Q' is preserved, the requirement that the two decay modes have 
different Q' charge implies that the two modes also have different B/L charge so that the set of 
particles comprising f\ and f 2 are different. This further implies that rescattering A\ 2 between 
the final states of the two decay modes involves both visible and hidden sector particles so that 
it is not possible to evade the bound (|66|) above. 

However, this bound may be evaded in theories where the global symmetry is of type S\ 
since then it is possible for fi and f 2 to contain the same set of visible sector particles. Rescat- 
tering then only involves the hidden sector and A\ 2 can be 0(1) in principle. There is still the 
requirement that the rescattering amplitude not wash-out the asymmetry once it is produced, 
but this is highly model dependent since the asymmetry may be produced at X" my , so that 
the washout is exponentially suppressed. Hence, in these theories e can be as large as 10 -2 for 
A n of Oil). 
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Figure 11: Cosmological phase diagrams similar to Figure [3 only depicting FO&D a and FI a along with other 
production mechanisms. Here the numerically evaluated contours correspond to £lh 2 arising from the sum of Y' 
and rj' . The left and right panels correspond to e = 10~ 8 and 10 -2 , respectively. 



Can these new mechanisms generate sufficient DM? This requires m'r]' = 4x10 10 GeV. For 
Asymmetric FI 



m'r}' m ~ 4 x KT 10 GeV 



10" 10 s\A 



12 



T 



a 



FI 



10- 6 J V 10- 3 I V 40 GeV 



m 



lOOGeV V sin0|A^ 2 | 
m J l^l 2 + \A 2 \ 2 



while for Asymmetric FO&D 



4 x 10" 10 GeV 



10- 8 (^)o|A 12 | 
(av) 10" 6 



FO 



10" 



m 



40 GeV ) 



100GeV\ smcplA^ 



m 



Since the last factor in these equations is always less than unity, Asymmetric FI requires a 
short lifetime r < 10~ 10 s |A 12 |/10 -6 and Asymmetric FO&D requires a small annihilation 
cross-section (av) < 10~ 8 (crt>) |A 12 |/10~ 6 . 

The symmetric yields Y' in (160]) and (|6T|) do not take re- annihilation into account. As has 
been discussed earlier in section I2.6[ if the (symmetric) yields for FO&D and FI are sufficiently 
large, then re-annihilation occurs giving rise to much smaller values for the final (symmetric) 
yields Vp 0&Dr and Vpj . This is crucial for the asymmetric yield to dominate the symmetric 
yield since e arises at the loop level and is expected to be small. DM will be dominated by the 
asymmetric component of X' only if (av)' is sufficiently large for re- annihilation to occur, hence 
reducing the symmetric component Y' while leaving the asymmetric component rj' unaltered. 
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Thus a DM asymmetry can be the dominant component of DM only in the re-annihilation regions 
for the symmetric abundance depicted in Figure [7J It is interesting that while FI dominates the 
DM density only in a narrow region of parameters as shown in Figure [7J asymmetric FI can 
dominate over a very wide region of parameters (assuming other conditions for generating the 
asymmetry are satisfied). 

4.2 Cosmological Phase Diagrams for DM Asymmetries 

Cosmological phase diagrams are shown in Figure [11] for non-zero values of e. The parameters 
held fixed are the same as in Figure [7J allowing a comparison of the asymmetric and symmetric 
DM production mechanisms. The contours are of Qh 2 , now arising from the sum of Y' and rj', 
and we have included purple and yellow regions denoting FO&D a and FI a phases, respectively. 
The left and right panels of Figure [TT1 correspond to e = 1CT 8 , 10~ 2 . Here the lower value of e 
arises from the non-thermalization bound in Eq. (|66|) and the upper value arises in theories of 
type 5*i where e is only bounded from above by the fact that in the perturbative regime, CP 
violating decays must always be at least a loop factor down from the CP respecting decays. 

In Figure [11] we see that the regions corresponding to FO&D and FI are essentially unchanged 
from Figure [7J while the regions for FO&D a and FI a are shown encroaching on the FO&D r and 
FI r regions. This is reasonable because these asymmetric DM production mechanism are only 
active when re-annihilation effects are maximal, so the symmetric component of the DM relic 
abundance is efficiently destroyed. As expected, for larger e, as depicted in the right panel of 
Figure [11] FI a becomes dominant at larger values of r, which is reasonable because in this case 
the CP phase in X decays is larger so the decay rate of X can be smaller while producing 
the same asymmetric DM abundance. On the other hand, at large e, FO&D a can dominate at 
smaller values of (o~v)'/ (cv), since the large CP phase in X decays allows for the abundance of X 
arising from FO to be smaller while producing the same DM asymmetry from late decays. There 
are no contours in the FO&D a region because this production mechanism gives an abundance 
independent of {cv)' and r. 

Note the behavior of the contours of Qh 2 depicted in Figure [TT] In particular, as soon as these 
curves enter the FI a region, they turn vertically upward. This is an indication that the abundance 
of DM contributed by Asymmetric FI is entirely independent of (av)'/(av), and instead only 
depends on r — thus, FI a behaves essentially identically to FI without re-annihilation effects. 

Next, let us consider the boundary between Asymmetric FI and Asymmetric FO&D. In 
particular, the abundance from FI a dominates the contribution from FO&D a when 1]' F1 > r]' FO&D , 
so 



r < a 



Mpi 




(67) 



FI 



m 2 Yl 



Cfo^fo 



fo&d 



m 
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Thus FO&D a dominates at large r, since then FI a is negligible, and for (av) <C (av ) , so there 
is a very large FO abundance of X. Notice that the right-hand side of Eq. (167|) depends on 
quantites which are fixed between the two panels of Figure [TTJ which is why the boundary 
between FI a and FO&D a is at the same value of r in both. 

FO&D a requires very small (crv), especially if e is small, and also dominates for large r where 
it might lead to signatures in big bang nucleosynthesis. On the other hand, FI a is independent 
of (av) and can successfully account for DM over a very wide range of r. The challenge for 
both asymmetric mechanisms is to construct theories where the CP violating phase in e can be 
measured in the laboratory. 

5 Conclusions 

The nature and origin of DM remains a deep mystery. A priori, there is no reason for DM 
to couple to the visible sector with interactions other than gravity, as only the gravitational 
properties of DM have been definitively verified through astrophysical probes. As such, one seeks 
theoretical motivations for why DM should be anticipated at other laboratory experiments, for 
example in direct detection, particle colliders, or indirect signals from cosmic rays. 

One popular theoretical justification is the so-called "WIMP miracle" , whereby DM is initially 
in thermal equilibrium with visible sector particles and undergoes FO, yielding approximately 
the correct DM relic abundance for weak scale masses and cross-sections. This setup has the 
advantage that DM generation is "IR dominated", and fixed wholly by the DM annihilation 
cross-section, which may be directly accessible at colliders. 

That said, it is clear that the single sector setup required for the WIMP miracle is only a 
very particular slice within the totality of theory space. In general, extra dimensional and string 
theoretic constructions motivate the existence of one or more hidden sectors which are very 
weakly coupled to our own visible sector. The inclusion of a weak-scale hidden sector, complete 
with its own set of particles, dynamics, and thermal history, substantially expands the allowed 
space of mechanisms for DM generation. Concretely, single sector FO is extended to two sector 
FO&D and FI, along with their re-annihilated and asymmetric cousins. 

Despite the proliferation of DM production modes, this broad two sector framework retains 
the IR dominated features lauded in the WIMP miracle. As we have shown, since both the visible 
and hidden sectors are initiated in a thermal state, the origin of DM and its final abundance 
are dictated entirely by the handful of quantities shown in Eq. (T5]). Remarkably, many of these 
quantities, namely m, m', (av), and r may be measured at particle colliders! This offers the 
exciting prospect that the origin of DM might be successfully reconstructed from collider physics 
in this enormous and theoretically motivated class of theories. We consider the details of such 
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an endeavor in a companion paper [8]. 
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